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ABSTRACT 

A  cylindrical  acoustic  resonator  which  has  been  developed  by  Dr. 
William  J.  Toulis  of  the  Signal  Propagation  Division  of  the  United  States 
Navy  Electronic  Laboratory  for  measurement  of  the  velocity  and  absorption 
characteristics  of  an  enclosed  fluid  medium  is  discussed.   The  cavity  is 
constructed  of  compliant  tubes,  thus  presenting  a  pressure  release  surface 
at  the  boundaries.   Theory  is  developed  indicating  that  the  resonant  fre- 
quency of  the  cavity  is  directly  proportional  to  the  velocity  of  sound 
in  the  medium  while  the  Q  of  the  resonator  can  be  used  to  determine 
attenuation. 
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1.   Introduction. 

The  success  of  sonic  echo-ranging  methods  depends,  in  part,  upon  the 
velocity  and  attenuation  profiles  in  the  transmitting  medium  between  the 
detecting  ship  and  the  target.   The  development  of  an  instrument  which 
could  directly  measure  these  properties  of  the  sea  would  be  of  great 
assistance. 

The  velocity  of  sound  in  sea  water  under  varying  physical  conditions 
has  been  investigated  both  experimentally  and  theoretically  by  many  per* 
sons.   For  the  experimental  determinations,  much  of  the  early  work  was 
pointed  toward  accurate  determination  of  depth  by  sonic  sounding  methods. 
Transmission  time  between  a  transmitting  and  receiving  hydrophone  was 
measured  in  an  area  where  salinity  and  temperature  were  as  uniform  as 
possible,  and  in  deep  enough  water  so  that  bottom  reflections  were  mini- 
mized.  In  this  type  of  determination,  measurement  of  the  distance  is 
the  limiting  factor. 

More  recently,  a  number  of  velocity  meters  has  been  developed  for 
this  purpose.   The  most  successful  for  laboratory  use  has  been  the  inter- 
ferometer, where  a  standing  wave  is  set  up  in  a  cell  by  utilizing  a 
crystal  and  reflector  or  two  reflectors.   This  then  provides  a  measure  of 
the  wave  length  and  hence,  the  velocity.   However,  no  completely  satis- 
factory operational-type  system  is  presently  available.!  The  most  widely 
used  device  to  determine  sound  velocity  profiles  is  the  bathythermograph, 
which  measures  temperatures  rather  than  velocities.   Its  use  depends  on 
the  fact  that  the  velocity  of  sound  is  largely  determined  by  the  tempera- 
ture of  the  sea,  and  only  to  a  lesser  extent  by  salinity  and  pressure. 
Urick^  obtained  good  general  agreement  between  velocity-depth  plots 
obtained  with  an  interferometer  and  those  computed  from  bathythermograph 
and  thermocouple  observations  for  sea  states  of  zero  to  three  and  depths 


J>H.  Sussman,  Study  and  Evaluation  of  ONR-Sponsored  Sound-Velocity 
Meters,  USL  Research  Report  No,  299,  1  March  1956 

2r.  J.  Urick,  An  Acoustic  Interferometer  For  the  Measurement  of 
Sound  Velocity  In  the  Ocean,  USNR  &  S  Lab.  Report  No.  S-18 


as  great  as  150  feet.   However,  the  bathythermograph  is  insensitive  to 
temperature  microstructures, J  and  as  such,  does  not  give  an  accurate 
indication  of  transmission  characteristics. 

A  great  deal  of  continuing  effort  has  been  given  to  the  theoretical 
determination  of  the  velocity  of  sound  from  the  physical  constants  of 
the  medium  as  given  by  the  theory  of  thermodynamics.  »"»'  Much  of  the 
work  is,  of  necessity,  empirical  in  nature  and  has  been  complicated  by  a 
lack  of  sufficiently  accurate  experimental  data. 8  However,  recent  im- 
provements in  laboratory  measuring  techniques  should  rapidly  improve  this 
situation. ^ 

Techniques  for  measuring  attenuation  of  sound  by  sea  water  are  even 
less  satisfactory.   The  absorption  at  low  frequencies  is  sufficiently 
small  to  make  measurement  difficult  without  considering  a  large  expanse 
of  sea.  Here  again,  determination  of  distance  is  a  limiting  factor. 
Further  isolation  of  the  specific  causes  of  attenuation  is  next  to  im- 
possible without  a  device  which  can  determine  the  absorption  character- 
istics of  small  samples. 

The  velocity  meter  under  discussion  depends,  for  its  operation,  on 
the  fact  that  the  resonant  frequency  of  a  water  column  in  a  particular 

3Urick,  ibid. 

^L.  Liebermann,  The  Effect  of  Temperature  Inhomogeneities  in  the 
Ocean  On  the  Propagation  of  Sound,  J.  Ac.  Soc.  Am.,  23,  pp  563  (1951) 

->S.  Kuwahara,  The  Velocity  of  Sound  in  Sea  Water  and  Calculation 
of  the  Velocity  for  Use  in  Sonic  Sounding,  Japanese  J.  Astron.  &  Geophys., 
Vol.  16,  No.  1  (1938) 

6n,  j.  Matthews,  Tables  of  the  Velocity  of  Sound  in  Pure  Water  and 
Sea  Water  for  Use  in  Echo-Sounding  &  Sound-Ranging,  Hydrographic  Dept. 
Admiralty,  London,  H.D.  282  (1927) 

?V.  A.  Del  Grosso,  The  Velocity  of  Sound  in  Sea  Water  at  Zero 
Depth,  NRL  Report  4002,  11  June  1952 

8r.  T.  Beyer,  Formulas  for  Sound  Velocity  in  Sea  Water,  J.  Mar. 
Research,  Vol.  13,  pp  113  (1954) 

9v.  A.  Del  Grosso,  E.  J.  Smura  and  P.  T.  Fougere,  Accuracy  of 
Ultrasonic  Interferometer  Velocity  Determinations,  NRL  Report  4439, 
6  December  1954 


mode  of  vibration  is  directly  proportional  to  the  velocity  of  sound  in 
the  water  column.   The  technique  employed  is  not  new  in  itself.  A  hard 
wall  cavity  resonator  was  developed  by  the  University  of  Michigan  in  1952.  " 
However,  this  device  was  extremely  sensitive  to  wall  cleanliness  due  to 
the  relatively  high  sonic  pressures  at  the  boundaries,  characteristic  of 
rigid  walled  devices.   Thus,  although  the  instrument  provides  an  accurate 
measure  of  the  velocity  when  used  in  the  laboratory,  it  is  very  difficult 
to  use  under  operational  conditions. 

Use  of  a  pressure  release  walled  cylinder  suggests  itself;  and,  in 
fact,  thin  walled  cylinders  have  been  used  for  measurement  of  the  atten- 
uation of  marine  sediments.  *» l   Recent  developments  in  the  field  of 
compliant  structures  by  Dr.  William  J.  Toulis  of  the  Signal  Propagation 
Division  of  the  United  States  Navy  Electronics  Laboratory,  San  Diego, 
California,  enabled  him  to  adapt  this  technique  to  the  present  operational 
device,  where,  by  measuring  the  resonant  frequency  and  the  "Q"  of  the 
cavity,  the  desired  characteristics  of  a  small  sample  of  water  can  be 
determined  to  a  high  degree  of  accuracy. 

The  purpose  of  this  paper  is  to  discuss  the  theoretical  aspects  of 
the  pressure  release  walled  cylindrical  cavity  constructed  by  Toulis, 
including  a  somewhat  systematic  study  of  the  loss  problem. 


*"R.  K.  Brown,  The  Development  of  an  Underwater  Sound  Velocity  Meter 
Using  a  Cylindrical  Resonator,  ProJ .  2021,  final  report,  30  April  1953, 
Eng.  Research  Inst.  Univ.  of  Michigan 

"W.  J.  Toulis,  Theory  of  a  Resonance  Method  to  Measure  the  Acoustic 
Properties  of  Sediments,  NEL  Report  676,  5  March  1956 

I^g.  Shumway,  A  Resonant  Chamber  Method  for  Sound  Velocity  and 
Attenuation  Measurement  in  Sediments,  Geophys.  Vol  XXI,  No.  2,  pp  305- 
319,  April  1956 
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2.  General  Wave  Equation. 

Before  proceeding  to  the  actual  analysis  of  the  velocity  meter,  it 
is  necessary  to  develop  the  equations  governing  the  transmission  of 
acoustic  waves  through  a  fluid  medium.   It  is  assumed  that  the  medium  is 
locally  homogeneous,  continuous  except  at  the  boundaries,  and  isotropic 
when  at  rest.   It  is  further  assumed  that  the  condition  of  the  medium  can 
be  fully  described  by  establishing  the  functional  dependence  of  pressure, 
vector  velocity,  density  and  temperature  on  time  and  the  space  coordin- 
ates. A  set  of  equations  consisting  of 

a)  the  continuity  equation  expressing  the  conservation  of  mass 

b)  a  force  equation  based  en  conservation  of  momentum 

c)  a  heat  exchange  equation  based  on  conservation  of  energy 

will  be  developed  in  their  complete  form  to  present  a  rigorous  point  of 
departure  for  making  the  necessary  approximations  for  the  study  of  small- 
signal  effects. 

To  derive  the  first  of  the  required  relationships,  consider  an 
elemental  volume  V  bounded  by  a  closed  surface  Sf.   The  amount  of  fluid 
crossing  the  surface  in  time  Atj  is  given  by  the  surface  Integral: 

outflow  m   At;  j  f  V*  A!S 

where  * 

V  ■  particle  velocity 
-f  a  point  fluid  density 

During  this  same  period,  the  quantity  of  fluid  originally  contained  in  V 
will  have  diminished  by  the  amount: 

loss  =  -At;  C  2-f  d  V 

where  the  negative  sign  is  selected  because  density  Is  a  decreasing 
function  of  time.   Since  mass  must  be  conserved,  net  outflow  must  equal 
the  loss  of  mass  and  the  two  expressions  can  be  equated.   The  surface 
integral  is  then  converted  by  use  of  the  divergence  theorem  giving: 

j>.(fv)  +f£]=0 

Since  the  integrand  is  continuous  and  the  volume  arbitrary,  this  may  be 
written  as: 

—•  +  v-  r-fvj  =  o  2.1 

which  is  one  form  of  the  well-known  equation  of  continuity. 
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To  develop  the  second  of  the  basic  relationships,  consider  Cauchy's 
14 


Considering  the  density  is  a  function  of  both  space  and  time,  the 
total  derivltive  of  ~P  with  respect  to  "D   can  be  written  in  vector 
notation  as: 

Of  =   ||  +  v-vf 

13       .-D 
utilising  the  notation  of  Stokes   where  W^T   denotes  the  material  de- 
rivltive and  is  defined  by 

Dt>     o  t   ^ 

Substituting  from  equation  2.1  into  2.2  gives 

— f  +•  fV*V  =  O  2.3 

Db 

which  is  Euler's  form  of  the  continuity  equation. 

To  develop  the 
first  law  of  motion: 

<m)  =  *  *  +  &? 

where: 

px  »  extraneous  force  vector,  a  function  of  position  and 

time 
Va.  =  velocity  component  of 
"Cu.*  stress  tensor. 
For  an  isotropic  medium  where  stress  is  a  linear  function  of  the  rate  of 
deformation,  the  stress  tensor  can  be  resolved  into  a  sum  of  a  scalar 
and  a  residual  stress  tensor  as  follows: 

-bA,  =  --P&Afr  +  UAfr         (^Af^t^A.)  2.5 

where 

^p  =  pressure 
IAa£»  residual  stress 
5a 4  «=  Kronecker  delta 
ad^a.l.=  hydrostatic  tension. 

13C.  True8dell,  The  Mechanical  Foundations  of  Elasticity  and  Fluid 
Dynamics,  J.  Rat.  Mech.  &Anal.,  1,  pp  125,  (1952) 

l^F.  V.  Hunt,  Notes  on  the  Exact  Equations  Governing  The  Propaga- 
tion of  Sound  in  Fluids,  J.  Acous.  Soc.  of  Am.,  27,  pp  1022,  Nov.  1955. 


To  a  first  approximation,  the  residual  stress  can  be  expressed  by 
the  linear  terms  of  a  power  expansion  in  the  viscosity  coefficients  as: 

U^=     XdKK8A£     +■    ZL/jdA^  2.6 

where: 

cU^  s  rate  of  deformation  tensor 
/J    =  dynamic  shear  viscosity  coefficient 
\  3  dilatatlonal  viscosity  coefficient 
and  dw   is  defined  by 

Combining  the  equations  2.5,  2.6,  and  2.7  and  substituting  into  2.4  yields 
the  vector  force  equation: 

-P(££)=  fF-f(7«7)V-Vf  +iX+2-yu)v(v-v) 

~/U7X(v^V)+  (v»-V)VX  i-  2-(V/>V) V     2.8 
+  V/JX(7X7). 
To  develop  the  third  of  the  fundamental  equations,  it  is  necessary 
to  establish  a  series  of  relationships  based  on  the  law  of  conservation 
of  energy.  The  total  energy,  b_T  ,  of  a  given  system  can  be  expressed  as 
the  sum  of  the  kinetic  energy,  p  .,  ,  and  the  internal  energy,  E.  ^ ,  where 

Ex  H  f -f  G  d  V  2*9 

'Xi   -  specific  volume  «  v-f 
£,  -  Internal  energy  density 
V  -  volume  of  the  system. 
The  material  rate  of  change  of  the  total  energy  is  then  the  sum  of 
the  rates  at  which  work  is  done  upon  the  system;  that  is: 

where 

"S"  «  heat  flux  vector. 
This  postulates  that  the  rate  of  doing  mechanical  work  can  be  divided 
into  two  portions:   that  arising  from  the  extraneous  force  and  that  aris- 
ing from  the  forces  of  material  continuity  expressed  in  terms  of  an 
equivalent  stress  vector  acting  on  the  surface  T>.   The  third  term  ex- 
presses the  rate  at  which  thermal  energy  enters  or  leaves  the  volume. 


Applying  the  divergence  theorem  end  utilizing  equetion  2.4  reduces  2.10 
to 

-?(§£)  =  ■&*♦<!*«  "Iff  '•» 

This  is  the  Fourier-Kirchhoff-C.  Neumann  energy  equation.  Using  the 
relationships  from  equations  2.5  and  2.7,  2.11  can  be  written  as 

Considering  the  internal  energy  density,  €-  ,  to  be  a  function  of  «S  end 
D  gives 

From  the  first  law  of  thermodynamics, 


2.13 


Thus,  equation  2.13  can  be  written  as 

-pfi_£c  -  -pT^  -  -^^-V  2.14 

Defining  ($>n  ,  a  viscous  dissipation  function  as 

and  eliminating  "P-=-S-  between  equations  2.12  and  2.14  gives:     # 

4>t25-    b    d>    -  3g  a  2.i5 

Next,  the  effects  of  radiation  and  conduction  must  be  considered. 
Since  the  divergence  of  "cj*   gives  the  energy  transferred  away  from  the 
volume  element,  it  must  account  for  energy  loss  by  either  conduction  or 
radiation. 

For  conduction,  the  Fourier  relation  Is  used: 

where  ^    / 

/j  •  thermal  conductivity. 
Taking  the  divergence  of  each  side  gives: 

XEa)cW   r/^v"T-(VT).(V/f)  2.16 

o  \a 

For  the  radiation  effects,  Newton's  law  of  cooling  states: 
»y™*   =fCv^Cr-Ti,)  2.17 


where 

Q  -  radiation  coefficient  introduced  by  Stokes,  Phil.  Mag. 

(4)  1,  pp  305-317  (1851) 

Cv  -  specific  heat  at  constant  volume  s  ^(AJL\ 

\*V  ~~  To)  -  local  temperature  excess. 

Before  proceeding,  the  term  T(^) appearing  in  equation  2.15  must 

be  expressed  in  terms  of  ~y"  »  \)     »    and  ""f  .   Following  the  methods  of 
.15 


+- 


2.18 


HuntAJ,  this  is  given  by 

T  2-  C?/cv 
O  -  coefficient  of  termal  expansion 

=  tey>T)f 

Combining  equation  2.18  with  2.15,  2.16  and  2.17  gives  the  third  of  the 
required  equations: 


-^vg?-^0 


VV-/<V  T 


-vT-v/r-i-fCv^Cr-  To)  -  4> 


n 


=    o 


2.19 


15Hunt,  op.  cit.,  pp  1023. 


3.   Small  Signal  Acoustic  Equations. 

The  small  signal  acoustic  equations  are  obtained  by  arbitrarily 
"linearising"  the  general  equations  of  motion.   For  this  purpose,  the 
dependent  variables  will  be  replaced  by  the  sum  of  their  equilibrium  or 
zero-order  values  and  the  first-order  variational  components.   Thus, 

f  =  "ft,  +  "ft  f  ■=    Jp0  4.  Jpt 

T  3  To  +  e,  v  =   o  +  v; 

On  the  assumption  that  -f^  ,  r~j~J)  and  too     are  constants  and  that  V 
is  itself  a  first-order  term  (i.e.~v0tr  O  )  then: 

vf  =  Vf,  v^=  Vf, 

vT=vOv  v-v  =  v-v, 

Substituting  these  relations  into  equations  2.1,  2.8  and  2.19  and  re- 
taining only  the  first-order  terms  gives: 

r~     +  -fD(v/N7.)  =0  3.1 


3.2 


3.3 


'I 

Po Cv9  ©i  ~  O  «> 

where  the  extraneous  force  has  been  omitted  since  it  primarily  affects 
the  equilibrium  configuration  and 

N  ~-  £L  "^*    /A'  )  tne  viscosity  number. 
Considering  JO   as  a  function  of  -P  and  ""J"1  gives: 


Since 


and 


(14?\     _     Co2- 


3.4 


equation  3.4  becomes 


0 


and  3.2  can  be  rewritten  as 

These  equations  (3.1,  3.3,  3.5)  take  Into  account  the  first-order  effects 
of  shear  and  dllatatlonal  viscosity  and  of  heat  conduction  and  their 
solution  Is  discussed  In  Section  5.   They  can  be  reduced  to  the  classical 
small-signal  wave  equation  by  neglecting  the  effects  of  heat  conduction 
and  viscosity  (/(=  Q  «/j  ■  \  «0)  and  assuming  adlabatlc  behaviour  of 
the  medium  giving 

f£  =#c»v-V 
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4.   Solution  of  the  Simple  Wave  Equation. 

Before  considering  the  more  complex  equations  derived  in  the  pre- 
ceeding  section,  a  solution  will  first  be  obtained  for  the  small-signal 
wave  equation.   The  results  obtained  will  serve  as  a  valuable  guide  when 
interpreting  the  solutions  obtained  in  later  sections. 

The  simple  wave  equation  is  given  by: 

^      J£  =  c  v   -p 

Where  V   is  the  Laplacian  operator  and  O  is  the  free  space  sound 

_  Afft, 

velocity.   Assuming  a  time  dependence  of  the  form  c     where  0"=-  CO 1-  AcX. . 
equation  4.1  takes  the  form  of  a  scalar  Helmholtz  equation: 

y*  +  Ka)-p  =-o  4.2 

where 

Using  the  method  of  separation  of  variables  with  the  Laplacian  operator 
in  cylindrical  coordinates  gives: 

J5."  +  J-  Jl'  +  ±,31'    +^+k^o       4.3 
I?       r  Ra    r^®     »  z 

where  p^u  denotes  4— B  ,  ^  denotes=r-^-a  ,  etc.   By  the  standard  argument 
for  the  method  of  separation  of  variables,  equation  4.3  may  be  rewritten 

"ja."  +  J__gl  +    '  -f-"  +  k*  =  -  -^  =  K/ 


0 


R    r  R  "^  t^  ® 

This  then  produces  two  equations: 

3-\  ±  £1  +  -J,  -S. 


^  -+  ^r-^~  -t^-M-   +  *   -   ^  4.4a 


Zu  +  K^Z  -    O 


4.4b 


where  equation  4.4b  is  the  familiar  differential  equation  of  simple 
harmonic  motion  with  the  solution: 

In  general,  Kj  can  be  considered  complex;  that  is 

Thus,  equation  4.5  may  be  rewritten  as 

2  =  (A^COS^E  +  B^^K^)  eT^         4.5a 
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Equation  4.4a  may  be  separated  utilizing  the  same  techniques  giving: 
©"-+  m^@  =  O  4.6a 

where  solutions  to  4.6  may  be  written  in  the  form: 

<®  =  Cm  cos(Vr,e  +  4-0  *•* 

Since  @  must  be  single  valued  in  0  ,YT\  must  be  an  integer. 

Equation  4.6b  is  recognized  as  a  form  of  Bessel's  equation  providing 
a  solution  (with  fn  as  an  integer)  of  the  form: 

where  JY^(yyV]is  a  Bessel  function  of  the  first  kind  and  NymKy^  is  a 
Bessel  function  of  the  second  kind16,  and  r\y-  =  V\  — K?  •   The  boundary 
conditions  which  must  be  satisfied  by  equations  4.5a,  4.7  and  4.8  are: 

a.  The  pressure  is  finite  throughout  the  interior  of  the  cylinder. 

b.  The  walls  are  assumed  to  be  pressure  release  surfaces  requiring 
that  Jjp  =  O  at  the  boundaries  of  the  cylinder. 

Selecting  a  coordinate  system  with  2  - O  at  the  bottom  of  the  cylinder, 
and  for  a  cylinder  of  radius  a  and  length  \__  ,  the  following  conditions 
apply: 

a.  Coefficient  of  the  cosine  term  in  equation  4.5a  must  be  zero 
since  ^p  is  zero  at  2.-  O* 

b.  Periodicity  in  Z.    is  given  by  K2  =  jt^"  where  #  »  1,2,3-  since  /p 
is  zero  at  H  =  L_  • 

c.  The  coefficient  of  the  second  term  of  equation  4.8  must  be  zero 
since  Bessel  functions  of  the  second  kind  cannot  be  present  with 
Jp     finite  at  the  origin.17 

d.  Since  the  selection  of  the  coordinate  system  is  purely  arbitrary 
the  0  =  0  axis  can  be  positioned  so  that  the  phase  angle  in 
equation  4.7  is  zero  (^«0)  with  no  loss  of  generality. 


^icLachlan,  Bessel  Functions  for  Engineers  -  2nd  Ed.,  Oxford  at 
the  Clarendon  Press  (1955) 

l'McLachlan,  op.  cit. 
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Thus,  the  solution  of  equation  4.2  takes  the  form: 


■c<o^ 


Acrt 


inCK^  cos(rne)X,(KvT)^    ^.9 


The  final  condition  to  be  satisfied  is  that  the  pressure  must  be  zero 
atr»(V   Thus 

Defining  a  quantity  ^>  vy\Y\  such  that 


4.10 


4.11 


then 

-  o 

where  P\  is  an  Integer  denoting  the  nth  zero  of  a  Bessel  function  of  order 
YT\   .  Further,  it  can  be  shown***  that  all  roots  of  equation  4.11  are 
real;  that  is,  Syr,n  must  be  real.   Thus 

is  real.   Since 


K\r  — 


4.12 


K 


NT 


-2.       2. 

K  -  K» 


then 


K^a  K-_(^y 


Substituting  the  complex  forms  for  Kj;  and  K  gives: 
Equating  the  real  and  imaginary  parts  of  equation  4.13  gives: 

Combining  equations  4.14a  and  4.14b  and  solving  for  K  gives: 


\2. 


If  b^,  is  assumed  negligable, 


K-^^+C^f 


4.13 


4.14a 


4.14b 


4.15 


4.16 


18 

A.  Gary  and  G.  B.  Mathews,  A  Treatise  on  Bessel  Functions,  Mac- 

millan  &  Co.  Ltd.  St.  Martin's  Street,  London,  pp  79,  1931. 
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where  the  negative  root  is  rejected  since  only  positive  frequency  is  being 
considered.  As  might  be  expected,  equation  4.16  is  the  equation  for  the 
wave  number  in  a  lossless  medium.   Substituting  for  K  and  \^,   in  equa- 
tion 4.16  gives: 


"Jsa.  =Wf +  PTf 


Therefore,  -f^Ynw »  tne  resonant  frequency  of  the  cavity  is 

«*»  =  iW  *  &S-Y ' 

It  ie  desirable  to  excite  a  mode  of  vibration  in  which  the  medium 
has  no  tangential  components  of  velocity.   For  a  medium  having  even  a 
slight  viscosity,  such  motion  at  the  walls  can  introduce  a  much  greater 
damping  of  free  vibrations  than  that  due  to  absorption  of  energy  by  the 
medium  itself.   Thus,  only  axially  symmetric  modes  will  be  considered 
(i.e.  m  *   0)  for  the  remainder  of  the  analysis.   For  the  lowest  axially 
symmetric  mode,  Sr^r^=  So|  »  x(0.7655)  and  ^  =  1.   This  corresponds  to 
a  pressure  antinode  at  the  midpoint  of  the  axis  of  the  resonator. 

If  attenuation  is  present,  its  effect  will  be  to  lower  the  resonant 
frequency  slightly  as  indicated  by  equation  4.15. 
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5.   Sound  Absorption  and  Dispersion  By  the  Medium. 

Investigation  of  the  general  mechanism  of  absorption  by  the  medium 
will  be  based  on  a  solution  of  the  first-order  equations  obtained  in 
section  3.   They  were 

H  +  -?0  (v«v)  ==  o  (31)     5>1 

^oCvdfjt^oCvCx-O/sIKv-v)  -  /(vze 

_        r  J  (3.5)   5.3 

t/^vx  (vxv)  =  o 

where  the  subscripts  indicating  first-order  effects  have  been  omitted. 

Differentiating  equation  5.3  with  respect  to  time  and  substituting 
from  5.1  f or  y^;  gives  ^ 

Considering  that,  according  to  Helmholtz's  theorem^  any  vector  can  be 
uniquely  separated  into  the  sum  of  a  gradient  of  a  scalar  and  the  curl  of 
a  divergence less  vector,  let 

V  =V(j)+Vx  A  5.5 

where 

CD  -  velocity  potential 
*/5^"  -  vector  potential 
V-7T-  0 
Substituting  equation  5.5  into  5.4  yields  the  following  scalar  and  vector 
wave  equations: 

*&  -4£»V*  +  ^0t.#  -/iN^(^)-0  5>6a 

-P0||    +/JVXVXA    =-0  5.6b 

It  is  interesting  to  note  that  equation  5.6b  which  governs  first-order 
vorticity  is  not  a  true  vector  wave  equation  but  a  vector  analogue  of  the 


l^Morse  and  Feshback,  Methods  of  Theoretical  Physics,  Vol.  1,  McGraw- 
Hill  Book  Co.,  pp.  53,  1953. 
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typical  homogeneous  diffusion  equation.   Since  no  source  terms  are  present 
in  5.6b,  this  implies  that,  if  first-order  vorticity  exists,  it  must  be 
generated  at  the  boundaries.   This  lack  of  source  terms  is  due  to  the 
technique  utilized  in  generating  the  first-order  equations.   Based  on  the 
assumption  that  the  equilibrium  values  of  the  state  variables  are  con- 
stant and  that  the  viscosity  coefficients  are  functions  of  these  state 
variables,  then  it  is  consistant  to  assume  that  V/Jq  «=  VAq  =  0  ■   There- 
fore, the  source  terms  involving  the  gradients  of  the  viscosity  coeffi- 
cients appear  as  second-order  rather  than  fits t-order  sources.   This  section 
deals  specifically  with  absorption  by  the  medium;  thus,  solution  of 
equation  5.6b  and  its  implications  will  be  reserved  for  section  6. 

Neglecting  the  effects  of  radiation  and  assuming  a  time  dependence 
\uSXs 
of  the  form  £     ,  equations  5.2  and  5.6a  become: 

V>P^\auVUN]    *    U>*-ff,(V-A60^/3-foe=.    O  5.7 


^24,[^oCv(^-i^]-^v^e  +  AoofGCv6  «  o 


5.8 


Solving  5.7  for  8  and  substituting  into  5.8,  and  introducing  the  short- 
hand notation:  * 
G 


and 


where 


gives 


-j — ^T»  a  frequency  dependent  term 


H 


GH 


,  which  is  non-frequency  sensitive 


— i   ii 

(|,<.^(j^AS(itVH)»-||2IF&CT*H-AH)JvO.O  5.9 

Equation  5.9  is  biquadratic  in  the  complex  propagation  constant  and  it  is 
convenient  to  rewrite  it  as 


(vV  K^)(V+-K,X)  4>  =  O 


5.10 


The  standard  quadratic  formula  can  be  used  to  determine  the  reciprocal 
squares  of  the  propagation  constants  as  -n 


|  +  AGr(l+TH)±£  +  AG(l-Ttf)j 
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For  the  frequency  range  under  consideration  with  water  as  the  medium,  G 
and  H  are  small  numbers  and  equation  5.11  can  be  expanded  by  power  series 
methods  to  yield: 


5.12a 


5.12b 


Neglecting  second-order  terms  gives: 

+■ 


KlSa  ±^U-^±{*l.  +(v-i\  A-l-^  ?  5.na 


13b 


Referring  again  to  equation  5.10,  we  see  that  the  solution  will  take  the 
form  of  two  waves  referred  to  as  types  I  and  II  whose  wave  numbers  are 
respectively  Ki   and  Kt_  •   Examination  of  equations  5.13  Indicates  that 
the  two  solutions  might  also  be  classified  as  "compressional"  and  "thermal" 
since  the  type  I  waves  represent  the  solution  to  the  ordinary  wave  equa- 
tion with  a  slight  damping  coefficient  and  the  type  II  waves  are  generally 
the  result  of  thermal  conduction. 

From  equation  5.13a  the  absorption  and  dispersion  measures  for  the 
type  I  waves  are : 


The  attenuation  thus  determined  is  the  Kirchhoff  approximation  for  low 
frequencies. 

The  wave  length  of  the  type  II  wave  (  \~2.TT(*%  ^J)  is  very  much 
shorter  than  that  of  the  type  I  wave.   Further,  the  type  II  waves  are 
very  rapidly  attenuated  at  the  frequencies  under  consideration  with  a 
decrement  of  "$—  — £>-  ,  and  the  velocity  of  propagation  through  the  medium 
is  given  by : 

~  'to  Cp" 
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6.   Wall  Losses. 

To  study  the  effects  of  wall  losses  on  the  system,  consider  the  force 
equation  derived  In  section  2  (2.9).   The  viscosity  coefficients  /u      and  \ 
are  assumed  constant,  and  effects  of  the  extraneous  force,  F,  are  assumed 
negligible.   Further,  since  only  small  amplitude  vibrations  are  to  be 
considered,  all  terms  involving  squares  of  the  velocity  components  are 
neglected.   Squat  ion  2.8  then  becomes 

-P(^)«  -V?   +(\^)7(7-\/)t/J7V  6.1 

Since  it  was  shown  in  section  5  that  the  thermal  waves  are  rapidly  atten- 
uated, the  effect  of  heat  conduction  at  the  boundary  will  be  neglected. 

Thus,  before  proceeding  further  in  the  analysis,  a  direct  relationship  be- 
on 
tween  pressure  and  velocity  must  be  established   .   Flow  of  mass  out  of 

any  volume  element  reduces  the  pressure  in  a  compressible  fluid.   In  an 

elastic  fluid,  if  compression  is  small,  the  rate  of  change  of  pressure  is 

proportional  to  the  divergence  of  the  velocity;  that  is 

-|£  -s  -"£>  VV  6.2 

B  -  bulk  modulus  of  elasticity 

Taking  the  time  derivitive  of  6.1  and  introducing  the  expression  given  in 

equation  6.2  for  -  £   gives: 

"ST/  

Again  defining  V  as 

-y^  v4>  +  VX/\  6.4 

and  substituting  equation  6.4  into  6.3  gives  the  two  equations: 

^Tz     *  ^7X7^  6-5a 

Assuming  a  time  dependence  of  the  form  6,     where  <T— W+ACX  ,  equations 
6.5a  and  6.5b  take  the  form: 


^^Morse  and  Feshbach,  Methods  of  Theoretical  Physics,  Vol.  1,  McGraw- 
Hill  Book  Co.,  pp.  162,  1953. 
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(v2fKl)^   =   0  '  6.6b 

where 


h 


2.  _   a -per 


k*  *  -±2± , 

Since  CD  satisfies  a  scalar  Uelmholz  equation,  the  form  of  solution  will 

be  the  same  as  equation  4.9,  that  is 

-r-  /  ,   \  Co5/       \     A  ®V 
<J>=SWK*2  J^(\<rr)  &^(vn  e)  1L  6.7 

where 


K* 

**»*  + 

Kr 

K* 

U 

A= 

l/2-i 

Jvr^ 

,(***} 

=  O 

Equation  6.6a  takes  the  general  form  of  a  vector  "wave"  equation.   Accord- 
ing 
by: 


21 
ing  to  Stratton   the  three  independent  vector  solutions  to  6.6a  are  given 


U=^V  6.8 

N    ^VXM  6.io 

where  0\    is  any  constant  vector  of  unit  length  and  y  is  a  scalar  func- 
tion satisfying  the  equation: 

From  6.8,  6.9  and  6.10, 

Since,  by  definition,  ^f- ft  =*  Q  »  6.9  and  6.10  must  be  the  complete  solution 
of  6.6a;  that  is,  * 

~K  ■=  "Ta   4-1T 

21J.  A.  Stratton,  Electromagnetic  Theory,  McGraw-Hill  Book  Co., 
pp  393,  1941. 
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cos   ^  6-14 


Equation  6.11  is  again  a  scalar  Helmholz  equation  and  thus  V   is 
given  by: 

Considering  only  the  zero  mode  (m  »  0) 

Substituting  equation  6.12  into  6.9  and  6.10  and  selecting  ~c\  -  O^  gives : 

N  =  h"(VKM) 

Taking  the  curl  of  6.13  and  6.14  gives 
From  the  symmetry  of  the  problem: 

vr(vya)  =   vr.(r,  l--i)  6.16a 

Thus,  the  rotational  component  of  the  vector  velocity  has  the  form 

r  s*D^g(>gcosy,  (&o]  6-17 

where  — • 

ba-    i1T/u         X--  1.2,;.--  - 

From  an  examination  of  equations  6.7  and  6.17,  it  can  be  seen  that  since 
Kv-  b  .  then  g  #  Ky  .   Thus     J0(g  a)  ^    O. 

The  irrotational  component  of  velocity  is  determined  by  taking  the 
gradient  of  6.7.   Thus 

20 


6.18 


A  rational  accounting  of  the  actual  wall  losses  is  well  beyond  the 
scope  of  this  discussion.   Derivation  of  even  a  reasonable  approximation 
for  the  pressure  field  existing  at  the  boundary  would  be  difficult  and 
tedious.   Even  if  these  equations  were  known,  the  mechanism  of  acoustical 
energy  transfer  to,  and  the  absorption  by  the  compliant  structure  is  not 
sufficiently  understood  to  make  application  feasible.   Thus,  it  is  nee* 
essary  to  assume  a  simple  model,  considering  that  the  solutions  obtained 
will  act  as  a  guide  in  interpreting  the  actual  behavior  of  the  resonator. 

As  a  first  approximation,  the  pressure  release  boundary  condition 
is  retained  with  the  additional  assumption  that  the  cylinder  is  bounded 
by  a  smooth  wall.   It  is  further  assumed  that,  for  a  viscous  medium,  there 
is  no  slipping  between  the  medium  and  the  wall.   This  implies  the  condi- 
tion that  the  tangential  components  of  velocity  will  be  zero  at  the 
boundary . 

Examination  of  equation  6.18  indicates  that  the  irrotational  com- 
ponents already  satisfy  this  condition,  i.e., 
Vr  *  O     oxt    2  =•  O  ,  "2  -  L 

v2  =  o       at      c  =  a  • 

Considering  the  rotational  components,  since  J,v§iCK-):£  0,  l/fcci  must  be  zero 
•o  that 

Further,  ^ b,oj  must  equal  zero  for  the  longitudinal  component,  Vs   ,  to 
satisfy  the  boundary  condition.   Thus,  for  this  simple  model,  there  is  no 
rotational  velocity  component,  and  no  first-order  vorticity  exists.   There* 
fore,  the  solution  obtained  in  section  4  is  unaltered. 

A  closer  approximation  to  the  actual  conditions  can  be  obtained  by 
assuming  that  pressure  approaches  some  small  finite  value  rather  than  zero 
at  the  boundary  (i.e. ,  J~0(  Kv^^fc  0  ,  bmVC^L  ~0  ).   In  this  case,  the 
coefficient  ~D 'o,g  in  equation  6.17  must  again  be  equal  to  zero.   However, 
from  the  two  conditions: 


two  equations  are  obtained 
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eP-KtSW^^-KKyV-V]    +    c[bg  s\v\bz  J"(Qr)j=o     6> 
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at     2r  ~  O     and    2  =  U.  }  and 


E[kzU)sV(^  Jo(«v^1  +    C^cosbz  7o{%o) 


=  0.  6.20 


Since  equation  6.19  must  hold  for  all  values  of  Vs  ,  it  follows  that 
K^b-  1-^.   Equation  6.20  can  then  be  rewritten  as 

y0(KroOf-  -^rC  3Usa) 

Equation  6.21  again  indicates  that  for  perfect  pressure  release  walls  (i.e., 

Xoi^v1^)  —  0  )  there  is  no  first-order  vorticity  and  no  viscous  loss  of 

energy  at  the  boundary. 

i 
As  the  boundary  deviates  from  the  ideal  case,  the  coefficient  O   in 

equation  6.21  increases  accordingly  and  vorticity  is  generated.   Thus, 

friction  arises  between  the  layer  of  fluid  adhering  to  the  walls  and  the 

immediately  adjoining  liquid  in  motion  causing  a  dissipation  of  energy  at 

the  walls. 

This  effect  can  be  evaluated  by  considering  that 

where  C-W\/n  is  a  small  complex  quantity  which  is  dependent  on  O   •   Re- 
ferring to  equation  4.13,  it  can  be  seen  that  the  addition  of  C  yy,^  will 
result  in  a  decrease  in  the  resonant  frequency  of  the  cavity  as  well  as 
an  increase  in  the  total  energy  dissipation. 

A  further  implication  of  this  analysis  is  that,  since  f^  4   O  at  T-  ^  , 
a  transmitted  as  well  as  a  reflected  wave  exists  at  the  boundary  and  some 
energy  will  be  lost  to  the  surrounding  medium  due  to  radiation. 
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7.  Measurement  of  Attenuation. 

Before  proceeding  to  apply  the  results  obtained  in  the  preceding  sec- 
tions to  the  actual  design  of  a  cavity,  some  method  of  measuring  the 
energy  absorption  by  the  enclosed  medium  must  be  determined.   The  quality 
factor  or  Q  of  a  resonant  cavity  may  be  defined  as: 

o  .  2jt  (energy  stored  in  cavity)  7  . 

"  energy  lost  per  cycle 

The  mechanical  energy  stored  in  an  incremental  volume  can  be  considered 

as  the  sura  of  the  kinetic  and  potential  energies;  that  is, 

ET  =   EK  +   EP  7.2 

where       E„  *±2ldv 

Therefore,  W,  the  total  energy  stored  in  the  cavity,  is  determined  by  tak- 
ing the  volume  integral  of  equation  7.2;  i.e., 

w«r^dv  i-  r  hpi^    j  y  7.3 

For  the  present,  wall  losses  are  neglected  and  all  attenuation  is  assumed 
to  occur  in  the  medium.   Selecting  a  time  dependence  of  the  form  z,        where 

0"  is  again  equal  to  ASXiCK)    and  substituting  into  equation  7.3  indi- 
cates  that  W  decays  at  a  rate  determined  by  £.       where  ^  is  the  atten- 
uation constant  in  nepers  per  second.   Thus,  the  energy  lost  during  one 
period  is  given  by: 

energy  lost  per  cycle  »  25  7.4 

Substituting  equations  7.3  and  7.4  into  7.1  gives 

n   airf  w  _  ..go  . 

or 

o(  s     --,     nepers/second  7.5 

To  convert  -X  from  a  time  to  a  space  dependence,  equation  7.5  is  multi- 
plied by  j/c  : 

^  g  '-^ZrC~  -  .r  ;-.   nepers/unit  length  7.6 

ci  C  C^     'Sl  w 

1 

where  K  is  defined  by  equation  4.15  as 


Y*^  (^^-^t^TA^*"] 


K  = 


Substituting  in  equation  7.6  and  salving  for  OC  gives: 

1 


-  Jl 


Y^  *  ( 


Thus,  with  the  assumption  that  all  energy  dissipation  occurs  in  the  medium, 
the  quality  factor  of  the  cavity  can  be  utilized  as  a  measure  of  the  at- 
tenuation constant  of  the  medium. 

In  the  case  of  an  actual  resonator,  the  effect  of  wall  loss  cannot 
be  neglected,  and  the  cavity  is  characterized  by  a  total  Q,  QT,  given  by 

_i_=    J_     +    _L_ 

where  Qwa  is  due  to  dissipation  by  the  medium  and  Qw  is  a  correction 
due  to  wall  loss.   Qm  ,  then,  is  the  value  which  must  be  utilized  in 
equation  7.7  to  measure  absorption  by  the  medium. 

Since  the  device  under  consideration  will  be  operating  at  relatively 
low  frequencies,  determination  of  the  correction  factor,  Q^  ,  is  greatly 
simplified. 

It  is  well  known that,  at  frequencies  below  10  k.c./sec,  attenuation 
of  sound  waves  in  fresh  water  is  considerably  less  than  attenuation  of 
sound  in  sea  water.   In  fact,  by  comparison,  attenuation  by  the  fresh 
water  can  be  neglected.   Thus,  if  the  cavity  is  immersed  in  distilled 
water  and  the  total  quality  factor  is  determined,  then 


Q-rd  Qwd 

to  a  close  approximation.  The  subscript  <A     is  used  to  indicate  measure- 
ment in  distilled  water. 

Q  may  also  be  defined  as 

Q-    -4r 

where 

«f  -  resonant  frequency 
Af   -  £>.  -  f , 

t^^j  -  frequencies  of  the  upper  and  lower  half  power  points. 
Therefore,  f  N         /   r  \ 


Considering  that  the  resonant  frequencies  in  fresh  and  salt  water  will  be 
only  slightly  different, 
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where 


Af. 


band  width  between  half  power  points  for  actual  cavity 
in  sea  water 


Ai-j  *  band  width  between  half-power  points  in  distilled  water 
X     -  resonant  frequency  in  sea  water. 


Substituting  this  value  into  equation  7.7  gives 

7.8 


tituting  this  value  into  equation  7.7  gives 


and  equation  7.8  can  be  used  to  experimentally  determine  attenuation  by 
the  medium. 
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8.   Cavity  Design. 

In  the  design  and  construction  of  a  sound -velocity  and  attenuation 
meter  based  on  the  principles  discussed  in  the  foregoing  analysis,  two 
general  problems  must  be  considered: 

a.  Design  of  a  cylindrical  cavity  presenting  essentially  a  pressure- 
release  boundary  over  the  frequency  range  under  consideration. 

b.  Development  of  the  "intelligence"  system  including  a  transducer 
to  excite  the  water  column  to  resonance,  a  receiving  transducer 
to  sample  the  generated  pressure  field  and  the  associated 
electronic  equipment  for  interpretation  and  display  of  the  de- 
sired information. 

As  a  first  step  in  the  design,  the  general  usage  of  the  data  obtained 
from  the  instrument  must  be  considered  to  determine  the  desired  resonant 
frequency  of  the  cavity.  Equation  4.17  can  then  be  utilized  to  calculate 
the  approximate  dimensions  of  the  cylinder.   Care  must  be  taken  to  select 
a  length  to  diameter  ratio  that  will  minimize  danger  of  mode  skipping. 

Next,  the  general  range  of  sound  velocities  to  be  measured  must  be 
specified.   This,  in  turn,  determines  how  broad  the  Q  of  the  cavity  must 
be  to  present  adequate  sensitivity  at  the  extremes  of  the  operating  range. 

With  these  parameters  in  hand,  the  formulas  developed  by  Toulis2^ 
can  be  utilized  to  determine  the  nature  and  spacing  of  the  compliant 
elements  for  construction  of  the  cavity  walls.   This,  then,  outlines  the 
steps  to  be  carried  out  in  the  solution  of  the  first  portion  of  the 
design  problem. 

In  selecting  and  mounting  the  transducers  for  exciting  the  water 
column  into  resonance  and  sampling  the  generated  pressure  field,  several 
factors  must  be  considered.  First  of  all,  since  the  material  generally 
used  for  enclosing  the  transducers  is  highly  absorptive,  it  is  desirable 
to  mount  them  in  a  relatively  low  intensity  sound  field.   Further,  it 
is  necessary  to  keep  separation  sufficient  to  minimize  direct  coupling. 

The  design  of  a  system  for  interpreting  and  presenting  the  data  ob- 
tained is  primarily  dependent  upon  the  actual  operating  requirements  and 
will  not  be  discussed. 

**W.  J.  Toulis,  Acoustic  Refraction  and  Scattering  With  Compliant 
Elements,  to  be  published. 
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9.  Conclusions. 

Examination  of  the  foregoing  analysis  indicates  that  the  solution 
to  the  simple  wave  equation  obtained  in  Section  4  describes,  to  a  very 
close  approximation,  the  conditions  existing  in  an  actual  resonator. 
From  Section  5,  it  can  be  seen  that  the  effects  of  heat  conduction  on 
the  resonant  frequency  and  attenuation  constant  are  negligible.   Section 
7  indicates  that,  with  walls  closely  approximating  an  ideal  pressure- 
release  surface,  the  effects  of  first-order  vorticity  can  also  be  neg- 
lected.  Thus,  for  design  of  a  system  employing  the  techniques  under 
discussion,  the  simplified  analysis  presented  in  Section  4  should  prove 
adequate . 

One  particularly  desirable  feature  of  this  system  is  that  it  mea- 
sures velocity  and  attenuation  simultaneously.   This  should  be  of  great 
assistance  in  determining  the  effect  of  air  bubbles  and  other  natural 
phenomena  on  sound  velocity. 

Another  desirable  feature  of  this  equipment  is  the  frequency  range 
employed.   Since  some  doubt  exists  as  to  the  validity  of  adapting  sound 
velocity  data  measured  at  megacycle  frequencies  to  much  lower-frequency 
applications,  an  accurate  device  operating  at  these  lower  frequencies 
should  be  a  welcome  addition. 
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APPENDIX  I 

CAVITY  CONSTRUCTION  AND  TEST 

A  cylindrical  acoustic  resonator  sound -velocity  meter  was  constructed 
by  Toulls  and  preliminary  tests  were  conducted  by  Toulls  and  the  author 
from  6  January  1957  to  15  March  1957.   This  appendix  outlines  the  methods 
of  construction  and  Instrumentation  utilized  and  describes  the  test 
results. 

The  pressure-release  walls  of  the  cavity  under  test  were  constructed 
by  deforming  commercially  available  aluminum  tubing  (Alcoa  3003-H-14,  1/2 
hard,  1"  O.D.,  0.036"  wall  thickness)  into  nearly  flat  or  elliptic  cross- 
section  and  welding  the  ends  to  prevent  flooding.   These  tubes  were  then 
secured  to  a  cylindrical  frame  with  2-1/2  inches  spacing  between  tubes 
forming  a  closed  cylinder. 

The  water  column  was  excited  into  resonance  by  a  cylindrical  barium 
tltanate  transducer  mounted  Just  below  the  compliant  tube  top  of  the 
cavity  near  the  cylinder  wall,  and  a  similar  transducer  utilized  to  detect 
the  pressure  field  was  mounted  directly  opposite. 

For  making  velocity  measurements,  the  output  of  the  receiving  trans- 
ducer was  amplified  and  applied  to  the  exciting  transducer.   A  self- 
excited  oscillating  loop  was  thus  formed.   The  loop  voltage  was  applied, 
through  a  frequency  multiplier  (X10),  to  a  Potter  150  KC.  Frequency 
Counter . 

For  attenuation  measurements,  the  amplitude  of  oscillation  is  an 
indirect  measure  of  Q.  Direct  determination  was  made  by  driving  the  cavity 
with  a  General  Radio  Beat  Frequency  Oscillator  (Type  1304-A)  and  measur- 
ing the  bandwidth  between  the  half -power  points. 

During  the  tests,  the  output  of  the  receiving  transducer  was  applied 
to  the  amplifier  through  a  calibrated  potentiometer.   The  amplifier  out- 
put was  monitored  and  the  gain  (potentiometer)  setting  required  to  main- 
tain constant  output  was  recorded  as  an  indication  of  the  attenuation.  A 
block  diagram  of  the  system  is  shown  in  Figure  1. 

To  evaluate  the  ability  of  the  instrument  to  measure  actual  sound 
velocity  at  sea,  a  number  of  tests  were  conducted  in  the  coastal  waters 
of  San  Diego,  California.  For  comparison  ^purposes  and,  particularly,  to 
isolate  the  effects  of  aallnity,  an  additional  set  of  runs  was  made  at 
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the  NEL  testing  facility  at  Sweetwater  Lake.   At  each  coastal  station, 
the  cavity  was  lowered  by  means  of  a  bathythermograph  winch,  and  the 
following  data  was  recorded  as  the  chamber  was  stopped  at  intervals  on 
the  ascent: 

a.  Resonant  frequency  as  indicated  by  the  Potter  counter. 

b.  Temperature  as  determined  by  a  thermistor  bead  attached  to  the 
side  of  the  cavity. 

c.  Gain  setting. 

d.  Depth  as  indicated  by  a  wire  metering  device. 

The  results  of  these  tests2-*  indicate  that  the  system  tested  pre- 
sents promise  for  development  as  a  survey  instrument  for  making  sound- 
velocity  profiles. 
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J.  E.  Wllle,  Preliminary  Test  of  A  Cylindrical  Acoustic  Resonator 
Sound-Velocity  Meter,  NEL  Tech.  Memo.,  to  be  published. 
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